Accepted in Astronomical Journal 

Preprint typeset using lAT^^ style cmulateapj v. 11/12/01 



O 

C/3 



VERY LOW-MASS STELLAR AND SUBSTELLAR COMPANIONS TO SOLAR-LIKE STARS 
FROM MARVELS I: A LOW MASS RATIO STELLAR COMPANION TO TYC 4110-01037-1 IN A 

79-DAY ORBIT 

.THAN D] 

Massimiliano Esposito'" '^, Scott W. Fleming^'^,B. Scott Gaudi", Luan Ghezzi^"-", Jonay I. 
Gonzalez Hernandez^*^, Brian L. Lee^, Keivan G. Stassun^'^^'^', Eric Agol\ Carlos 
Allende Prieto'^'*, Rory Barnes^ Dmitry Bizyaev^**, Phillip Cargile^, Liang Chang^, 
Luiz N. Da Costa^" ", G.F. Porto De Mello^"'^"', Bruno Femeni'a'''*, Leticia D. 
Ferreira^"'^'', Bruce Gary'^ Leslie Hebb'^ Jon Holtzman^", Jian Liu^, Bo Ma^, Claude E. 
Mack IIL', Suvrath Mahadevan'^'^', Marcio A.G. Maia^"'", Duy Cuong Nguyen^, Ricardo 
L.C. Ogando"' ", Daniel J. Oravetz", Martin Paegert^, Kaike Pan", Joshua Pepper'', 
Rafael Rebolo Basilic Santiago"' "\ Donald P. Schneider'*'^^, Alaina C Shelden", 



(N 
O 

. Audrey Simmons", Benjamin M. Tofflemire^-^", Xiaoke Wan^, Ji Wang^, Bo Zhao 



Accepted in Astronomical Journal 

ABSTRACT 



fSJ . TYC 4110-01037-1 has a low-mass stellar companion, whose small mass ratio and short orbital period 

are atypical amongst solar-like (Tg// < 6000 K) binary systems. Our analysis of TYC 4110-01037-1 
reveals it to be a moderately aged (<5 Gyr) solar-like star having a mass of 1.07 ±0.08 Mq and radius 
of 0.99 ±0.18 R0. We analyze 32 radial velocity measurements from the SDSS-III MARVELS survey 

00- 

as well as 6 supporting radial velocity measurements from the SARG spectrograph on the 3.6m TNG 
(~| , telescope obtained over a period of ^2 years. The best Keplerian orbital fit parameters were found to 

have a period of 78.994 ±0.012 days, an eccentricity of 0.1095±0.0023, and a semi-amplitude of 4199±11 
m s~^. We determine the minimum companion mass (if sini = 1) to be 97.7 ±5.8 Mj„p. The system's 
companion to host star mass ratio, >0.087 ±0.003, places it at the lowest end of observed values for 
short period stellar companions to solar- like {Tf.ff < 6000 K) stars. One possible way to create such a 
I system would be if a triple-component stellar multiple broke up into a short period, low q binary during 

the cluster dispersal phase of its lifetime. A candidate tertiary body has been identified in the system via 
single-epoch, high contrast imagery. If this object is confirmed to be co-moving, we estimate it would be 
. , a dM4 star. We present these results in the context of our larger-scale effort to constrain the statistics 

^ • of low mass stellar and brown dwarf companions to FGK-type stars via the MARVELS survey. 

\0 ■ Subject headings: stars: individual (TYC 4110-01037-1) 

■ 1. INTRODUCTION influence our understanding of the process of planetary 

C^l ■ „, r , J- J. -u r 1 4_ 11 1 and star formation (Burgasser et al. 2007: Stamatellos 

OObservations ot the mass distribution ot substellar and „ T,ri ■; .i ^n^r, . i n , i . i 

^ . , , r . • r- 1 1 , • & Whitworth 2009: Kraus et al. 2011; Sahlmann et al. 

fvj low- mass stars as a function of orbital separation, us- r,^-,-,\ i i r i i / ■ i 

. -^riu- -JI 2011). Early results from the Kepler transit search pro- 

T-H ■ mg a variety of techniques, provide key constraints which .,. , ., 

... i r J gram, which is sensitive to short to moderate period com- 
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panions, show the size distribution of candidate compan- 
ions increases towards smaller planets, reaching a maxi- 
mum at a few R0 (Borucki et al. 201 la, b). Numerous large 
radial velocity surveys are constraining the prevalence of 
Jovian-mass giant planets at close and intermediate orbital 
separations (e.g. California & Carnegie teams, Marcy & 
Butler 2000; AAPS, Tinney et al. 2001; CORALIE, Udry 
et al. 2000; HARPS, Mayor et al. 2004) and the prevalence 
of brown dwarf (BD) and very low mass stellar compan- 
ions at close and intermediate separations (Marcy & Butler 
2000; Mayor et al. 2001; Vogt et al. 2002; Patel et al. 2007; 
Lee et al. 2011). At much larger orbital separations, high 
contrast imaging surveys inform our understanding of the 
mass distribution of planets (Kalas et al. 2008; Marois et 
al. 2008, 2010; Lagrange et al. 2010) and BDs (Thalmann 
et al. 2009; Biller et al. 2010; Janson et al. 2011; Wah- 
haj et al. 2011), while the mass distribution of very low 
mass stellar companions has been explored by high con- 
trast imaging (Kraus et al. 2011), spectroscopic (Duquen- 
noy & Mayor 1991), and interferometric (Raghavan et al. 
2010) studies. 

Many of the results from these surveys have found in- 
teresting, and at times conflicting, trends motivating sub- 
sequent programs to investigate their origin. In the BD 
regime, a deficit of short period (a < 5 AU) companions 
to solar-type primary stars has long been designated the 
apparent "brown dwarf desert" (Marcy & Butler 2000). 
Investigations of the BD frequency at larger orbital sep- 
arations (e.g. Gizis et al. 2001; Metchev & Hillenbrand 
2004; McCarthy & Zuckerman 2004; Metchev & Hillen- 
brand 2009; Kraus et al. 2011) have led to conflicting as- 
sertions as to whether the BD desert extends to larger 
orbital separations. Multiplicity studies of low-mass stel- 
lar companions to solar-like stars have found evidence of 
a unimodal period distribution with peak periods rang- 
ing from ^180 years (Duquennoy & Mayor 1991) to '--^300 
years (Raghavan et al. 2010). Interestingly, recent work by 
Metchev & Hillenbrand (2009) has presented tentative ev- 
idence that the companion mass function of brown dwarf 
and low-mass stellar companions around solar-like stars 
could be represented by a universal function. 

At the short-period (P <100 days) tail of solar-like 
{Teff < 6000 K) multiplicity investigations, many stud- 
ies have reported a paucity of confirmed low-mass stellar 
or brown dwarf companions with mass ratios {q = Mc/M^,) 
<0.2 (see e.g. Pont et al. 2005; Bouchy et al. 2011), lead- 
ing Raghavan et al. (2010) to suggest that short-period 
companions to solar-like stars prefer higher mass ratios, 
although there is disagreement in the literature regarding 
this trend (see e.g. Halbwachs et al. 2003). Burgasser et 
al. (2007) relate that the short-period brown dwarf desert 
seems to extend into the M dwarf regime. Raghavan et 
al. (2010) also note that the majority of the companions 
surrounding solar-like stars with periods <100 days are 
triple systems, and thus potentially indicative that such 
systems experienced orbital migration (Bate et al. 2002). 
Short-period low mass ratio binaries are much more com- 
monly observed around slightly more massive F-type stars 
(Bouchy et al. 2005; Pont et al. 2006; Beatty et al. 2007; 
Bouchy et al. 2011). Bouchy et al. (2011) have suggested 
that low q companions can form at or migrate to short 
orbital periods around a wide range of stellar primary 



masses, but suggest that they may not survive around G 
dwarfs and lesser mass primary stars. 

The Multi-object APO Radial Velocity Exoplanet 
Large-area Survey (MARVELS), one of the three sur- 
veys being executed during the Sloan Digital Sky Survey 
(SDSS) III (Eisenstein et al. 2011), is a four-year program 
which is monitoring the radial velocities of ^3,300 V=7.6- 
12 FGK type dwarfs and subgiants. As described in Lee et 
al. (2011), the target selection strategy attempts to impose 
minimal and well-understood biases on targets' ages and 
metallicities; hence, the survey provides an ideal, statisti- 
cally robust means to explore the mass distribution of sub- 
stellar and very low mass star companions over orbital pe- 
riods of <2 years from a relatively low-biased target sam- 
ple. In anticipation of a statistical analysis of global trends 
in the population of BDs and low mass binary companions 
identified by the survey, we are performing detailed and 
careful characterization of the fundamental parameters of 
these companions and their host stars (see e.g. Fleming et 
al. 2012 discussion of MARVELS-2b). The first paper in 
this series was Lee et al. (2011), which presented an analy- 
sis of a short period brown dwarf surrounding the F9 star 
TYC 1240-00945-1. Analogous detailed characterization 
of individual systems has been shown to be particularly 
important in refining the radii of Kepler candidate planet- 
candidates to be more Earth-like than assumed (Muirhead 
et al. 2011). The advantages of meta-analyzing such well- 
characterized systems is also demonstrated in Bouchy et 
al. (2011), who were able to begin to explore the mass- 
radius relationship from the planetary to BD to very low 
mass star regimes as well as the mass ratio of companions 
as a function of primary mass. 

In this paper, we present a detailed analysis of the fun- 
damental properties of the solar-like star TYC 4110-01037- 
1 (hereafter TYC 4110) and report on the discovery of a 
very low mass stellar companion associated with the sys- 
tem. In Section 2, we describe the spectroscopic and pho- 
tometric data which were used for this analysis. We de- 
termine accurate fundamental stellar parameters for the 
star in Section 3 and describe the basic properties of the 
very low mass stellar companion in Section 4. Finally, we 
discuss the implications of these results in Section 5. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. SDSS-III Radial Velocity Data 

Our primary radial velocity (RV) observations of TYC 
4110 were obtained during the first two years of the SDSS- 
III MARVELS survey, which uses a dispersed fixed-delay 
interferometer (Ge et al. 2009) on the SDSS 2.5m telescope 
(Gunn et al. 2006). A total of thiry-two observations were 
obtained over the course of ~2 years. Each 50 minute 
observation yielded two fringing spectra (aka. "beams") 
from the interferometer spanning the wavelength regime 
~500-570 nm with R ~12,000. Lee et al. (2011) describe 
the basic data reduction and analysis leading to the pro- 
duction of differential RVs for each beam of the interfer- 
ometer. We combined these beams for each observation 
set using a weighted mean. As described in Fleming et al. 
(2010), we scaled up the RV errors by a "quality factor" 
(Q — 5.67 for TYC 4110) based on the rms errors of the 
other stars observed on the same SDSS-III plate as TYC 
4110. 
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A summary of the relative amplitude RV measurements 
obtained for TYC 4110 with MARVELS is presented in 
Table 1. 

2.2. TNG Follow-up Radial Velocity Data 

Supporting RV observations were obtained with the 
3.6m Telescopio Nazionale Galileo (TNG) using its SARG 
spectrograph (Gratton ct al. 2001). The O'.'S x 5'.'3 slit pro- 
vided R ^57,000 spectroscopy between 462-792 nm. Wc 
obtained six spectra with an iodine cell (IC), to provide 
high precision radial velocities (Table 1), and one with- 
out the IC to serve as a stellar template. The data were 
reduced using standard IRAF routines, and RVs were mea- 
sured using the IC technique (Marcy & Butler 2000) . Each 
of 21 SARG spectral orders between 504-611 nm were di- 
vided in 10 pieces, and RV calculations were derived from 
each of the 210 resulting pieces. Based on a gooclness- 
of-fit indicator, the best 158 (75%) pieces were selected. 
Following a 2-a clip, the remaining RV measurements were 
combined with a weighted average to produce the RV mea- 
surements quoted in Table 1. 

2.3. APO Spectroscopic Data 

Two R~31,500 optical (-3,600 -10,000A) spectra of 
TYC 4110 were obtained on UT 2010 September 29 (HJD 
2455468) with the Apache Point Observatory 3.5m tele- 
scope and ARC Echelle Spectrograph (ARCES; Wang et 
al. 2003), to enable accurate characterization of stellar fun- 
damental parameters. The two spectra were obtained us- 
ing the default 1'.'6 x 3'.'2 slit and an exposure time of 1200 
seconds. A ThAr lamp exposure was obtained between 
these integrations to facilitate accurate wavelength cali- 
bration. The data were processed using standard IRAF 
techniques. Following heliocentric velocity corrections, 
each order was continuum normalized, and the resultant 
continuum normalized data from each observation was av- 
eraged. The final spectrum yielded a signal-to-noise ratio 
(S/N) of ~175 at ~6500 A. 

2.4. HAO Photometric Data 

We obtained absolute photometry of TYC 4110 using 
the Hereford Arizona Observatory (HAO), a private facil- 
ity in Southern Arizona (observatory code G95 in the lAU 
Minor Planet Center). HAO employs a 14-inch Meade 
LX200GPS telescope equipped with a SBIG ST-IOXME 
CCD. Observations in B, V, and Ic filters were made on 
2011 January 15 and 2011 February 10. A total of 22 Lan- 
dolt standard stars and 9 secondary standards based on 
Landolt star fields SA98 and SA114 (Landolt 1992) were 
observed at several airmass values similar to the airmass 
for TYC 4110. The target's magnitude was calculated by: 

Mf = Mfo - 2.5 logio(i^//5) - {K'f ■ m) + {Sf ■ C) (1) 
where M/ is the magnitude for each filter /, M/o is a con- 
stant determined from the standard stars, Ff is star flux 

using a large photometry aperture, g is exposure time, 
is zenith extinction coefficient for each filter, m is air mass, 
Sf is the star color sensitivity (determined from the stan- 
dard stars) , and C is star color (B- V) . Solutions for (B-V) 
were obtained by iterating V and B magnitudes from ini- 
tial values 2-3 times. The resultant absolute photometry 
for TYC 4110 is summarized in Table 2. 



2.5. SuperWASP Photometric Data 

We analyzed 5 epochs of broadband optical photometry 
(400-700 nm) of TYC 4110, obtained between 2006 April 
7 and 2008 April 14, from the SuperWASP public archive 
(Butters et al. 2010). Aperture photometry of the 160 in- 
dividual observations available in the archive, each taken 
with a 30 second integration time, was computed via the 
SuperWASP pipeline. Further details about the observa- 
tional design and data reduction pipeline for SuperWASP 
can be found in Pollacco et al. (2006). We find the Super- 
WASP photometry (exhibits no statistically significant ev- 
idence of variability (error-weighted RMS ~0.675%) over 
the time-scales sampled by these data. For example, a lin- 
ear fit to the entire dataset yields negligible variation in 
flux with a best-fit slope of 0.040 ± 0.030% day-^ We 
also detect no evidence of a transit. We do caution how- 
ever that the sampling of these data is sparse. 

3. TYC 4110-01037-1: the star 

3.1. Fundamental Stellar Properties 

We analyzed moderate resolution spectroscopic data 
from the ARCES spectrograph using two separate anal- 
ysis techniques to extract fundamental stellar parameters 
for TYC 4110. We refer to these different pipeline re- 
sults as the "lAC" (Instituto de Astrofisica de Canarias) 
and "BPG" (Brazilian Participation Group) results, as de- 
scribed below. 

3.1.1. "lAC" Analysis 

We derive equivalent widths (EWs) of Fe I and Fe II 
lines with the code ARES (Sousa et al. 2007), using 
an initial linelist with 263 Fe I and 36 Fe II lines given in 
Sousa et al. (2008), and use the rules in Sousa et al. (2008) 
to modify the rejt parameter in ARES according to the 
signal-to-noise ratio of each spectrum. In addition, we set 
the ARES parameters smoother = 4, space = 3 and mini- 
line — 2. The parameter lineresol was modified according 
to the resolving power of each spectrum. 

The stellar atmospheric parameters were computed us- 
ing the code StePar (Tabernero et al. 2012). This code 
employs the 2002 version of the MOOG code (Sneden 
1973), and a grid of Kurucz ATLAS9 plane-parallel model 
atmospheres (Kurucz 1993). StePar, iterates until the 
slopes of A (Fe I) vs. x and A (Fc I) vs. log (EW A^^) 
are equal to zero, while imposing the ionization equilib- 
rium condition A (Fe I) = A (Fe II). A 2-(t rejection of the 
EWs of Fe I and Fc II lines is performed after a first deter- 
mination of the stellar parameters, and then the StePar 
program is re-run without the rejected lines (see Tabernero 
et al. 2012, for further details). 

For the ARCES spectrum of TYC 4110, 198 Fe I fines 
and 24 Fe II lines remain after clipping. These are 
used to derive T^s = 5879 ± 25K, log{g) = 4.53 ± 0.18, 
[Fe/H] = -0.02 ±0.05, and Vmicro = 0.932 ± 0.038km 8"^ 

Internal uncertainties were also derived for each stel- 
lar parameter. The uncertainty of Vmicro was obtained 
by varying this parameter until the slope of the linear re- 
gression of A (Fe I) versus log (EW )^~^) was equal to its 
standard deviation. The uncertainty of Teff was deter- 
mined by changing this parameter until the slope of the 
linear regression of A (Fe I) versus % was equal to its stan- 
dard deviation. The uncertainty of Vmicro was also taken 
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into account when calculating the uncertainty of Tcff . The 
uncertainty of log (g) was obtained by varying this pa- 
rameter until the difference between the mean abundances 
from Fe I and Fe II were equal to the standard deviation 
of the latter. The contributions from Tgff and Vmicro were 
included. Finally, the uncertainty of [Fe/H] is a combina- 
tion of the standard deviation of the Fe I abundance and 
the variations caused by the errors in Tgff , log {g) and ^t, 
all added in quadrature. 

3.1.2. "BPG" Analysis 

We assume LTE and use the 2002 version of MOOG 
(Sneden 1973), along with the ID plane-parallel model at- 
mospheres interpolated from the ODFNEW grid of AT- 
LAS9 models (Castelh & Kurucz 2004). Initially, a list of 
^ISO isolated and moderately strong (i.e., 5 < EW < 120 
mA) Fe I and Fe II lines was compiled using the Solar Flux 
Atlas (Kurucz et al. 1984), the Utrecht spectral line compi- 
lation (Moore et al. 1966), and a Ganymede ARCES spec- 
trum with S/N = 400. The values for the central wave- 
lengths and line excitation potentials were taken from the 
Vienna Atomic Line Database (VALD; Kupka et al. 1999). 
We also multiplied the van der Waals damping parameter 
"C6" by a factor of two, following Holwcger et al. (1991). 

The EWs of these lines were automatically measured 
in the solar spectrum from fits of Gaussian profiles using 
the task hplot in IRAF. The quality of the measurements 
was checked by performing two tests. First, since the line 
depth is expected to be a linear function of the reduced 
EW (EWA~^) for non-saturated lines, we eliminated lines 
that did not follow a linear relation, using a 2-<j clipping. 
The second test is based on the fact that the shapes of 
the lines are essentially determined by the instrumental 
profile at the APO resolution (~31,500). Since the reso- 
lution is approximately constant over the entire spectrum, 
we expect the quantity FWHMA"^ to be approximately 
constant for lines of the same species. We therefore per- 
form a linear fit to this relation and eliminate lines that 
exhibit 2-cr deviations. 

After these tests, the final solar line list contained 91 
Fe I and 11 Fe II lines. Solar gf values were derived 
for all these lines using a solar model atmosphere with 
the following parameters: Tefi = 5777 K, log {g) = 4.44, 
[Fe/H] = 0.0 and Vmicro = 1-0 km s~^. The adopted solar 
abundance for iron is A(Fe) = 7.50 (Asplund et al. 2009). 

EWs for 102 Fe lines were measured in the TYC 4110 
APO spectrum and checked with the tests described above. 
Tgff and Vmicro were iterated until zero slopes were found 
in the plots of A(Fe I) versus x and log (EW A"'^) , respec- 
tively; i.e., until the individual Fe I line abundances were 
independent of excitation potential and reduced EWs. The 
surface gravity was iterated until A(Fe I) = A(Fe II), i.e., 
until the same average abundances were given by Fe I and 
Fe II lines. At the end of this iterative process, a consistent 
set of atmospheric parameters (Tos, \og{g), [Fe/H], and 
Vmicro) was obtained for the star. Note that the metallicity 
is simply given by [Fe/H] = A (Fe) — 7.50, where 7.50 is the 
solar iron abundance taken from Asplund et al. (2009). At 
this point, any lines with abundances that deviated more 
than 2-cr from the average were removed and the above 
iteration was repeated until convergence was achieved. 

We derive T^s = 5878 ± 49 K, log{g) = 4.43 ± 0.17, 



[Fe/H] = 0.00 ± 0.06 and Vmtcro = 1-00 ± 0.08 km s^^ 
based on the ARCES spectrum. The final line list after 
rejections contained 60 Fe I and 8 Fe II lines. The internal 
uncertainties are calculated in the same was as the "lAC" 
analysis above. 

3.1.3. Final Stellar Parameters 

We determined the mean values for Toff, log (g), [Fe/H] , 
and Vmicro by combining the results from the lAC and 
BPG analyses via a mean, weighted by the inverse of the 
internal variances. For each parameter, we add in quadra- 
ture a systematic error of 18 K, 0.08, 0.03 and 0.02 km s~^ 
for Teff, log(.g), [Fe/H], and Vmicro, respectively. These 
systematic errors are calculated based on the weighted 
standard deviation of the weighted means of each parame- 
ter using 18 spectra of 13 stars (7 MARVELS targets and 
6 stars with well-known atmospheric parameters). These 
stars span Teff from 5200-6500 K, log (g) from 4.0-4.7, 
[Fe/H] from -0.5 - -1-0.5 and Vmicro from 0.3 - 1.8 km s~^. 
The final stellar parameters are Teff = 5879 ±29 K, logg 
= 4.48 ±0.15, [Fe/H] = -0.01 ±0.05, and Vmicro = 1-40 
±0.05 km s-i (Table 2). 

To check these parameters using a different observa- 
tional technique, we constructed a spectral energy distri- 
bution (SED) for TYC 4110, using the near UV to 2MASS 
(Cutri et al. 2003) and WISE (Wright et al. 2010) IR pho- 
tometry compiled in Table 2 (see Figure 1). These data 
were fit with a NextGen model atmosphere (Hauschildt et 
al. 1999), and we limited the maximum line-of-sight extinc- 
tion to be Ay to < 0.28 from analysis of dust maps from 
Schlegel et al. (1998). The resultant parameters, Teff = 
6000 ±200 K, log5 = 4.0 ±1.0, [Fe/H] = 0.0 ±0.3, and Ay 
= 0.20±0.08, all agree to within 1-a of the results found 
via analysis of our moderate resolution spectroscopy. We 
also fit these SED data by constraining the Tef / and log g 
values to those derived from our spectroscopic analysis. 
The resultant fit (reduced = 2.62; see panel b Figure 1) 
provides a more robust estimate of Ay, 0.16 ±0.04. Us- 
ing this total extinction estimate, and adopting a BCy of 
-0.19 ±0.02 (Cox 2000), we estimate the distance to TYC 
4110 to be 125.1 ±4.6 pc (Table 2). 

3.2. Stellar Mass and Radius 

Using the spectroscopically determined values of Teff, 
log 5, and [Fe/H] (Table 2), we determined the mass and 
radius of TYC 4110 using the empirical Torres et al. (2010) 
relationship. We find M* = 1.07 ±0.08 Mq and R* = 0.99 
±0.18 R0 (Table 2). Our quoted errors include contri- 
butions from the uncertainties in our fundamental stel- 
lar parameters, as well as the scatter in the Torres et al. 
(2010) relationship {aiogm — 0.027 and aiogr = 0.014) and 
correlations of the best-fit coefficients from Torres et al. 
(2010) added in quadrature. We did not include covari- 
ances between Teff, logg, and [Fe/H] in this error analy- 
sis; however, our final quoted uncertainties for these values 
do include a systematic error term that conservatively en- 
capsulates any covariance between these parameters. We 
did however perform a Monte Carlo simulation of our spec- 
troscopically determined stellar parameters and the Torres 
et al. (2010) relations, and found a stellar mass and radius 
consistent with the aforementioned values. 
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3.3. Evolutionary State 

We assess the evolutionary state of TYC 4110 by com- 
paring its spcctroscopically measured fundamental stellar 
parameters against a Yonsei-Yale stellar evolutionary track 
(Demarque et al. 2004) for a = 1.07 M© star having 
[Fe/H] = -0.01. This is done in Figure 2, where the shaded 
region depicts deviations in the evolutionary track which 
would be expected from a l-o" (0.08 Mq) change in the 
assumed stellar mass, while circles denote different time 
stamps in the track. TYC 4110 lies near a predicted age 
of <5 Gyr in Figure 2; we therefore conclude that it is a 
main sequence dwarf star. The lack of any detectable Ca 
II H and K emission in our ARCES spectra ( \og{Rjjj^') 
~ -5.1) indicates the star is relatively inactive, and thus 
qualitatively consistent with the evolutionary state we de- 
rive. 

3.4. Systemic and Rotational Velocity 

We computed the absolute radial velocity of TYC 
4110 by cross-correlating the six epochs of SARG spec- 
tra against a Solar spectrum. After removing the radial 
velocity contribution at each epoch induced by the pres- 
ence of TYC 4110's companion (Table 1), we determine 
the systemic velocity of TYC 4110, Vsystemic-, to be 25.9 
±0.2 km s-i (Table 2). We also compared these SARG 
data to broadened versions of Kurucz ATLAS synthetic 
spectra to constrain the rotational velocity, Vrotsini, of 
TYC 4110. After considering a range of macroturbulence 
values, we find Vrotsini <3 km s^^, which is slightly below 
the level of instrumental broadening present in these data 
(-5.3 km s"^). 

4. TYC 4110-01037-1'S COMPANION 

4.1. Binary Companion Detection 

Analysis of post-pipeline processed MARVELS data in- 
volves searching for periodic behavior which is consistent 
with Keplerian orbital motion, and filtering out "contam- 
inant" RV signals which do not arise from the presence 
of a companion. One of the first steps in this process 
is to compute a Lomb-Scargle (LS) periodogram (Lomb 
1976; Scargle 1982) and search for periodic signals which 
significantly exceed a conservative false alarm probability 
assessment. There have been several implementations of 
this method (Press & Rybicki 1989); we followed Gumming 
(2004) which uses the comparison of the values between 
a sinusoidal fit at a given frequency and a fit to the mean 
to generate the power term (see e.g. Cumming's equation 
2). This form was chosen to be fully extensible by adding 
linear terms, harmonic terms, or even a fully Keplerian 
fit. We stepped through frequency space using steps and 
a search window as appropriate for our data sampling, as 
described in Press et al. (1992). In order to interpret the 
significance of this power spectrum, false alarm probabili- 
ties (FAP) were calculated using the techniques in Baluev 
(2008). 

Data from the current MARVELS pipeline have system- 
atic effects present which can mimic a companion signal. 
We can significantly mitigate this issue by taking advan- 
tage of the fact that the MARVELS instrument observes 
60 stars at a time. As a result, we can search for peri- 
odic signals from all stars on a given plate, and use any 



detected signals to characterize systematics in our data. 
This is done by taking the sum of the power for each fre- 
quency across the entire plate. We then remove the high- 
est power from each frequency (so that actual companions 
do not skew the average), excluding any power associated 
with the candidate companion, and compute the average. 

Figure 3 indicates TYC 4110 exhibits a strong ~79-day 
period signal. The companion to TYC 4110, hereafter re- 
ferred to as MARVELS-3B, does not match any strong 
periodic signature identified in the power spectrum of all 
stars located on the same SDSS-III plate (bottom panel; 
Figure 3), and is therefore not caused by any known sys- 
tematic artifact in our instrument or reduction pipeline. 

4.2. Radial Velocity Fits 

Radial velocities derived from MARVELS and SARG 
data were fit with the EXOFAST code (Eastman et al. 
2012) to extract detailed Keplerian orbital parameters. We 
first performed an independent fit of the MARVELS data, 
and re-scaled the MARVELS error bars to force the prob- 
ability of P(x^)) = 0.5. Since we did not have enough 
SARG data points to perform an independent fit solely 
on these data, we then fit the combined (MARVELS + 
SARG) data, and re-scaled the SARG errors to force P(x^) 
= 0.5. The fitting of the combined dataset and re-scaling 
of the SARG errors were iterated until a convergent so- 
lution was achieved. The MARVELS and (MARVELS -|- 
SARG) data yielded consistent results to within l-cr of the 
fit errors. We also computed the of the SARG data 
about the MARVELS-only fit, and the resultant value (x^ 
= 8.47 for 5 degrees of freedom) indicates such a fit would 
only happen by chance —13% of the time. We hereafter 
only consider the combined MARVELS + SARG fit. 

The raw MARVELS and SARG radial velocities were 
computed on independent, relative scales. We determined 
that the best offsets to simultaneously analyze these data 
about a zero-point of m s^^ were 1345±24 m s^^ (SARG) 
and 8583 ±12 m s^^ (MARVELS); the RVs quoted in Ta- 
ble 1 have had these offsets applied. The RV errors listed 
in Table 1 include the re-scaling factors described above. 
The inclusion of an additional linear term was explored in 
the fitting process, but there is no compelling evidence for 
an acceleration due to another companion, with a best-fit 
linear slope of 0.0961^; m s ^ day ^ observed, corre- 
sponding to a 2.4-cr deviation. 

As seen in the full RV curve (Figure 4) and in the phase- 
folded curve (Figure 5), these data are described by a 
78.994 ±0.012 day period, a moderately elliptical (e = 
0.1095 ±0.0023) orbit, and a semi-amplitude K = 4199 
±11. Full fit parameters for TYC 4110 are presented in 
Table 3. To search for evidence that MARVELS-3B was 
a transiting system, we computed a LS periodogram of 
the available SuperWASP photometry, but found no evi- 
dence of variability at the ~79 day period of MARVELS- 
3B above a level of ~0.8%. 

We determined the mass for MARVELS-3B using: 

(MeSinz)3 _ j^3(x_g2)(3/2)p 

^1^1^ = (5.953 ± 0.047). 10-Me (3) 

and the fit parameters compiled in Table 3. This yields 
a minimum companion mass (if sini = 1) of 97.7 ±5.8 
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Mjtip, which places MARVELS-3B slightly above the gen- 
erally accepted brown dwarf upper mass limit of 80 Mj„p 
and into the low-mass star regime. The minimum mass 
ratio, q, of the companion to the primary is 0.087 ±0.003. 

4.3. Binary Companion Mass 

The true mass of the companion depends on the incli- 
nation i of the orbit, which is unknown. However, we can 
estimate the posterior probability distribution of the true 
companion mass, given an isotropic distribution of orbits, 
and adopting a prior for the distribution of the compan- 
ion mass ratios. We proceed to do this using a Monte 
Carlo method, following the methodology described in de- 
tail in Fleming et al. (2010) and Lee et al. (2011), which 
we briefly summarize here. We combine the posterior dis- 
tribution of orbital parameters e, and P obtained from 
the Markov Chain Monte Carlo (MCMC) fit to the ra- 
dial velocity data, with an estimate of the joint distribu- 
tion of the primary mass and radius obtained using the 
spectroscopically-dctermined T^/f, logg and [Fe/H] com- 
bined with the Torres et al. (2010) relations, accounting 
for all sources of uncertainty in the measured values and 
the relations themselves. We draw values of cosi from a 
uniform distribution. The values of K, e, and P determine 
the mass function (M^ sini)'^/(M, -I- Mc)^, and then the 
value of i along with the primary mass determines Mg. 
Finally, we appropriately weight the resulting distribution 
of Mc by our prior on the mass ratio q. 

As described in Section 1, the mass ratio of companions 
around G dwarfs is not well constrained by current obser- 
vations and MARVELS-3B likely lies in a relatively under- 
populated region of mass ratio parameter space. Neverthe- 
less, we consider several different priors on the companion 
mass ratio which we suggest arc reasonable given current 
observations (see e.g. Grether & Lineweaver 2006), of the 
form: dN/dq oc q+^, dN/dq oc q~^, and dN/dq = constant. 
We note that massive companions are ruled out by the 
lack of a statistically significant infrared excess in the SED 
(Figure 1) and the lack of a secondary component in our 
high-resolution optical spectra (see discussion below). We 
include this constraint by weighting the resulting distri- 
bution of Mc by exp[—0.5{AK/AK„iax)'^], where AK^ax 
is the upper limit on the excess flux (in magnitudes) in 
the K band, and AK is the excess flux contributed by a 
companion of mass Mc and a primary of mass M* , as de- 
termined using the Baraffe et al. (1998) solar-metallicity, 
Y = 0.275, 1 Gyr mass-magnitude relations (note that we 
could have adopted any isochrone in the range of 1-10 Gyr 
with negligible difference) . We caution the reader that this 
specific constraint does not fully and uniformly represent 
the prior for every possible configuration of the compan- 
ion, but rather is a reasonable, simplying set of conditions. 
We did not observe a statistically significant IR excess flux 
in TYC 4110's SED (Figure 1); thus, we used the 3-f7 stan- 
dard deviation (0.06 magnitudes) of the Ks-band data for 
AKmax- We note that the IR fiux contribution from the 
wider separation candidate tertiary companion, discussed 
in more detail in Section 4.4, is less than this AKmax and 
therefore does not influence our mass estimate. 

The resultant cumulative distributions of the true mass 
are shown in Figure 6, and we summarize the median mass 



for each of these priors, along with the resultant mass ratio, 
in Table 4. Since these probabilistic median values rely on 
the simplifying assumptions we have made, we have not as- 
signed formal confidence limits to these values, to prevent 
over-interpretation of their robustness. For each prior, we 
note that the resultant median mass is in the M dwarf re- 
gion with a low (<0.2) mass ratio, and companions more 
massive than 0.5 Mq are ruled out at the 95% probability 
level. 

Finally, we also analyzed the ARCES spectrum of TYC 
4110, obtained at a phase of ~0.71, to search for evidence 
of MARVELS-3B. We first note that we detected no evi- 
dence that the system was a SB2, which suggests q < 0.65 
(Halbwachs et al. 2003), i.e. MARVELS-3B is less massive 
than a mid K-type dwarf. Figure 7 illustrates the red opti- 
cal difference spectrum computed by subtracting the spec- 
trum of a G dwarf with very similar fundamental stellar 
parameters, HD 153458, from TYC 4110. Both the known 
binary companion to TYC 4110 (MARVELS-3B) and the 
candidate tertiary companion (sec Section 4.4) could fill in 
the spectral features identified with red dots in Figure 7, 
leading to positive deviations in the difference spectrum, 
although differentiating the precise amount of any contri- 
bution from the secondary versus the candidate tertiary is 
not possible. While calibrating our usage of this technique, 
we found that minor mismatches between the properties 
of the reference and science spectra could produce notica- 
ble subtraction residuals in line wings (as seen in Figure 
7), which suggests that only strong, repeatable deviations 
in these difference spectra should be interpreted as real 
contributions from companion(s). Using spectra from the 
Pickles (1998) library, we estimate that a M3-type com- 
panion would have contributed a ^2% flux enhancement 
to the system. The difference spectrum does not exhibit 
>3-cr deviations above the level of the Poisson noise (0.8%) 
and continuum normalization uncertainties present in the 
data, which sets the lower mass limit for MARVELS-3B 
that can be ascertained from these specific data. 

4.4. Candidate Tertiary Companion 

Triple star systems are a relatively common outcome 
of the star formation process (see e.g., Tokovinin 2004). 
To further assess the multiplicity of TYC 4110, we ac- 
quired adaptive optics (AO) images on 2012 January 7 
using NIRC2 (PI: Keith Matthews) at the Keck Obser- 
vatory in natural guide star mode. Our initial data set 
consisted of nine dithered images taken in the K filter. 
Inspection of the raw frames showed evidence for a faint 
candidate companion located to the south-west of TYC 
4110. Figure 8 shows the fully processed K ' image. 

We measured an accurate position for the candidate 
companion using the technique described in Crcpp et al. 
(2012). We first fit Gaussian functions to the stellar and 
companion point-spread-functions to locate their centroids 
in each frame. The primary star was not saturated in any 
of our dithered images. We then correct for distortion in 
the NIRC2 focal plane (narrow camera mode) using the 
publicly available software provided by the Keck Observa- 
tory astrometry support page The results are averaged 
and the uncertainty in the separation and position angle is 
taken as the standard deviation, taking into account un- 
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certainty in the plate scale and orientation of the array by 
propagating these errors to the final calculated position. 
Adopting a plate scale of 9.963 ± 0.006 mas pixel"^ and 
instrument orientation relative to the sky of 0.13° ± 0.02°, 
as measured by Ghcz et al. (2008), we find a companion 
separation and position angle of p = 986 ± 4 mas and 
PA = 218.1° ± 0.3° respectively. 

Upon noticing the companion, we obtained additional 
images in the J and H-bands to facilitate characteriza- 
tion. Our aperture photometry indicates that the object 
has red colors: A J = 4.219 ± 0.104, AH = 3.940 ± 0.032, 
and AK' = 3.805 ± 0.027 mags. Table 5 lists its appar- 
ent magnitude as measured relative to the primary star, 
taking into account the combined light from each source. 

Both the (J-H) = 0.55 and (J-K') = 0.75 colors indi- 
cate a spectral type of ~M3V (Leggett et al. 2002). As- 
suming the candidate is situated at the same distance 
as the primary, we find that the absolute magnitudes, 
Mj = 8.10 ± 0.21, Mh = 7.55 ± 0.32, Mk = 7.36 ± 0.13, 
are each consistent with a ~ O.25M0 star when compared 
to the Girardi et al. (2002) evolutionary models, which 
corresponds to a main-sequence spectral type of ~dM4. 
The relations from Table 5 of Kraus & Hillenbrand (2007) 
yield the same result (dM4) , though with a possibly lower 
mass estimate of ~ O.2OM0. 

Given the RA and DEC and distance of TYC 4110, the 
a priori likelihood of detecting a background star within 
I'.'O is -0.8%. With a proper motion of [-0.20, -99.40] 
mas yr~^, a time baseline of several months will be suffi- 
cient to assess its association with the primary star. We 
note that the detection and mass constraints placed on 
MARVELS-3B in this paper hold regardless of whether or 
not this candidate tertiary companion is confirmed to be 
co-moving with the primary star. 

5. DISCUSSION 

We discuss some of the derived properties for 
MARVELS-3B in the context of previous studies of low- 
mass companions to solar-like stars. Both Duquennoy & 
Mayor (1991) and Raghavanet al. (2010) demonstrate that 
companions to solar-like stars having orbital periods <12 
days are circular; however, companions having orbital peri- 
ods similar to that of MARVELS-3B, ~79 days, exhibit ec- 
centricities from 0-0.6 (Duquennoy & Mayor 1991; Mayor 
ct al. 1992; Raghavan et al. 2010). MARVELS-3B's mod- 
est eccentricity of ~0.11 therefore is clearly consistent with 
that observed for other similar period companions. 

Arguably the most distinctive feature of TYC 4110 is its 
extremely low mass ratio, q, of >0.087 ±0.003, given its 
relatively short —79-day period. As illustrated in Figure 
9, previous statistical investigations of binarity in solar- 
like {Tgff < 6000 K) stars have found evidence that the 
short period brown dwarf desert extends in mass toward 
the low mass star regime (Burgasser et al. 2007; Bouchy et 
al. 2011; Sahlmann et al. 2011). A similar desert of short- 
period, low q companions to low-mass K and M dwarf stars 
has also been observed (Figure 9; see also Burgasser et al. 
2007). The studies of companions aroimd solar-like stars 
by both Duquennoy & Mayor (1991) and Raghavan et al. 
(2010) report no firm detection of low q binary companions 
with orbital periods <100 days. The closest short period 
analog observed by Raghavan et al. (2010) (see e.g. their 



Figure 17) is a multiple (> 2 components) q —0.23 system 
with an orbital period of —50 days, while the closest short 
period binary is a —1-day period object with a g of —0.4. 
All other low q binaries in their sample have orbital peri- 
ods >5000 days. Mayor et al. (2001) also report a sparse 
number of low q binary companions to solar-like stars, but 
do not quantify the orbital periods of these objects. 

The absolute "dryness" of the low mass ratio, short pe- 
riod desert for solar-like (Te// ^ 6000 K) illustrated in Fig- 
ure 9 is a subject of active investigation and debate in the 
literature. OGLE-TR-122b (Pont et al. 2005), Kepler-16b 
(Doyle et al. 2011), and vB 69 (Bender & Simon 2008) are 
short period, low q, solar-like binary systems that populate 
this parameter space. The statistical frequency of objects 
like OGLE-TR-122b and Kepler-16b have not been quan- 
tified by the OGLE or Kepler surveys; hence, it is plau- 
sible that they are indeed rare. We note that there have 
been suggestions the mass ratios of binary companions sur- 
rounding F7 to K-type stars are more broadly distrubuted 
and that no clear low mass ratio desert exists (see Figure 
8 in Halbwachs et al. 2003). However, we have included 
all of the companions surrounding G-type (0.8 M0 < M < 
1.1 Mq) and K- and M-type primaries tabulated in Halb- 
wachs et al. (2003) (but omitting data from their survey 
which was not tabulated in print or online) in Figure 9, 
and still clearly see a short period, low mass ratio compan- 
ion deficit. Clearly, additional observational studies which 
precisely characterize the stellar properties of the host star 
and characterize the properties of their binary compan- 
ions are needed to perform an accurate meta-analysis of 
the dependence of companion mass ratios as a function of 
orbital period and host star mass. These studies should 
make careful note of their observational biases to facili- 
tate a more robust cross-correlation of results compiled 
from different surveys. A future meta-analysis would also 
particularly benefit from having current and future studies 
fully tabulating the fundamental properties of the binary 
systems they investigate. 

As noted in Section 1, short-period, low q companions 
have been reported around stars slightly more massive 
that the Sun (F-type stars; Bouchy et al. 2005; Pont et 
al. 2006; Beatty et al. 2007; Bouchy et al. 2011)! Bouchy 
et al. (2011) propose that short-period low q companions 
might form around a wide mass range of stars, including G- 
type stars, but suggest that weaker magnetic disk braking 
during the early formation history of F-type stars might 
transfer less angular momentum to their companion bod- 
ies, thereby preventing catastrophic decays of their orbits. 
Conversely, Bouchy et al. (2011) propose that a stronger 
disk braking in young G-type stars might distribute more 
angular momentum to their companion bodies, caiising 
(initially) short-period companions to migrate inwards and 
become engulfed by the primary. While intriguing, this 
proposed evolutionary scenario would clearly benefit by a 
more robust assessment of how "dry" the short-period, low 
q desert is around GKM dwarfs as compared to F dwarfs. 

The precise mass ratio of MARVELS-3B is not known 
due to the unknown inclination of the system, thus our 
observations only set a lower limit of q of >0.087 ±0.003. 
However, as demonstrated in Section 4, our Bayesian anal- 
ysis of the system using 4 plausible prior assumptions all 
indicate the likely median mass of MARVELS-3B is a dM 
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star, yielding mass ratios of 0.087 < q < 0.149 (Figure 6 
and Table 4). Like OGLE-TR-122b (Pont et al. 2005) and 
Kepler-16b (Doyle et al. 2011), MARVELS-3B therefore 
seems likely to be an outlier to the mass ratio-period rela- 
tionships commonly observed for both solar-like stars and 
lower- mass dM stars (Figure 9). 

In this context, we note that analyses of the period dis- 
tribution of exoplanets have revealed a deficit of such bod- 
ies having orbital periods of 10-100 days, aka the "period 
valley" (Udry et al. 2003; Jones et al. 2003). Wittcnmycr 
et al. (2010) suggest that this observed deficit is real for 
the giant planet population (M > 100 M0), but that any 
deficit of lower mass planets (10-100 M^) in this period 
regime might be the result of selection effects. One pos- 
sible explanation for the dearth of giant planets in 10-100 
day orbits is a decrease in the amount of orbital migra- 
tion which such objects experience (see e.g. Trilling et al. 
1998). 

Migration could also play a role in setting the observed 
period distribution of the low-mass stellar binary regime. 
The orbital evolution of binaries has been explored com- 
putationally, and stellar accretion, the interaction between 
binaries with their natal gas disks, and interactions be- 
tween triple components can influence these systems (see 
e.g. Bate et al. 2002; Stamatellos & Whitworth 2009; 
Kraus et al. 2011). Bate et al. (2002), for example, found 
such processes were successful at producing short-period 
binaries in high mass ratio systems (q > 0.3), but it is 
uncertain as to whether these specific simulations could 
produce short-period, low q systems like MARVELS-3B. 

N-body simulations of the early evolution of stellar clus- 
ters which include instantaneous gas removal (Moeckel & 
Bate 2010) arc beginning to better reproduce the number 
of unequal mass solar-like binaries observed (Duquennoy & 
Mayor 1991; Raghavan et al. 2010). Potentially relevant 
for MARVELS-3B, Moeckel & Bate (2010) showed that 
one simulated triple system comprised of a 0.96 and 0.73 
M0, 2 AU separation binary with a 0.21 M0 tertiary at 12 
AU, broke up into a tight (0.1 AU separation) binary sys- 
tem comprised of the 0.21 and 0.73 Mq components (e.g. 
a mass ratio, q, of 0.29). We therefore speculate that it is 
possible that the short period, low q MARVELS-3B binary 
was initially part of a tertiary (or larger) system with much 
different initial orbital parameters, and only achieved its 
final orbital configuration following the dispersal of the 
cluster in which it formed. Although our analysis of avail- 
able SuperWASP photometry and MARVELS-hSARG ra- 
dial velocity data exhibited no evidence of the presence of 
additional, long-period (> 2 years) bodies in the system, 
our single-epoch detection of a candidate tertiary body in 
the system via high contrast imaging could support this 
interpretation. Additional epochs of imagery should be 
pursued to establish if this body is co-moving with the 
TYC 4110 system. Although beyond the scope of this 
paper, future simulations of binary migration and orbital 
evolution during cluster dispersal should explore the limits 
and frequency at which they can reproduce short period, 
low mass ratio binaries for solar-like stars. 

6. CONCLUSIONS 

We present a detailed analysis of the fundamental prop- 
erties of the solar-like star TYC 4110-01037-1 and its very 



low mass stellar companion. This analysis was performed 
in the context of our long-term goal of performing a de- 
tailed statistical analysis of global trends in the population 
of well vetted and characterized BDs and low mass binary 
companions identified by the MARVELS survey. We find: 

• TYC 4110-01037-1 is a <5 Cyr solar-like star char- 
acterized by Teff = 5879 ±29 K, \ogg = 4.48 ±0.15, 
and [Fe/H] = -0.01 ±0.05. We determine the stel- 
lar mass to be 1.07 ±0.08 Mq and stellar radius to 
be 0.99 ±0.18 R©. 

• MARVELS-3B is a >97.7 ±5.8 Mj^p (Msini) com- 
panion to TYC 4110-01037-1, which follows a mod- 
erately elliptical (e = 0.1095 ±0.0023), 78.994 
±0.012 day orbital period. 

• The mass ratio, q, of the companion to the primary 
is > 0.087 ±0.003. MARVELS-3B therefore resides 
in a short period, low q desert analogous to the short 
period brown dwarf desert. 

• We speculate that MARVELS-3B might have ini- 
tially formed in a tertiary system with much differ- 
ent orbital parameters, and achieved its present day 
configuration following the dispersal of the cluster 
in which it formed. A candidate tertiary body has 
been identified via single-epoch, high contrast im- 
agery. If this object is confirmed to be co-moving, 
we estimate it would be a dM4 star. 
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Table 1 

Summary of Observed Radial Velocities 



HJD 


Instrument 


RV 








km 


km 


2454811.815473 


M 


-2.869 


0.053 


2454812.936774 


M 


-2.386 


0.049 


2454816.915562 


M 


-0.841 


0.060 


2454840.865816 


M 


3.859 


0.046 


2454842.900506 


M 


3.614 


0.093 


2454843.814123 


M 


3.496 


0.057 


2454844.776592 


M 


3.342 


0.045 


2454845.779836 


M 


3.153 


0.049 


2454867.724621 


M 


-2.229 


0.072 


2454868.791455 


M 


-2.527 


0.048 


2454869.785190 


M 


-2.786 


0.070 


2454901.671359 


M 


1.393 


0.068 


2455105.974063 


M 


-2.454 


0.046 


2455135.969729 


M 


0.417 


0.039 


2455141.825627 


M 


2.508 


0.045 


2455142.892315 


M 


2.744 


0.060 


2455143.904775 


M 


2.978 


0.058 


2455144.904270 


M 


3.187 


0.042 


2455161.844226 


M 


3.115 


0.047 


2455199.808399 


M 


-4.240 


0.059 


2455201.786093 


M 


-4.002 


0.040 


2455202.835898 


M 


-3.854 


0.042 


2455280.657591 


M 


-4.064 


0.046 


2455287.644203 


M 


-2.150 


0.052 


2455463.908559 


M 


3.816 


0.073 


2455470.923509 


M 


4.025 


0.042 


2455472.000492 


M 


3.942 


0.047 


2455498.003488 


M 


-1.714 


0.066 


2455500.992667 


M 


-2.521 


0.059 


2455516.486551 


S 


-4.209 


0.061 


2455516.579785 


S 


-4.118 


0.033 


2455542.731058 


M 


3.798 


0.049 


2455545.783002 


M 


4.117 


0.053 


2455553.524971 


S 


3.677 


0.036 


2455556.930997 


M 


3.175 


0.076 


2455580.495813 


S 


-2.702 


0.034 


2455666.464340 


S 


-4.051 


0.037 


2455698.384891 


S 


3.375 


0.041 



Note. — A summary of relative radial velocities obtained with the MARVELS (M) and SARG (S) spectrographs. The quoted aRv errors 
for the MARVELS data were first uniformly scaled by a "quality factor" Q =5.67 (see Fleming et al. 2010), based on the rms of the other stars 
observed on the same SDSS-III plate. As described in Section 4, the errors of the MARVELS and SARG data were then re-scaled to force our 
RV fits to have P(x^) = 0.5, in an iterative maimer. Note that zero point ofTsets have been applied to the RVs compiled in column 3 to force 
the RVs to vary about m s~^. 
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Table 2 

Stellar Properties of TYC 4110-01037-1 



Parameter 


Value 


Uncertainty 


Note 


a (2000) 


06 54 12.3 




Hog et aL (1998) 


<5 (2000) 


+60 21 31.4 




Hog et al. (1998) 


NUV 


15.715 mag 


0.016 


GALEX (Morrisscy et aL 2007) 


B 


11.159 mag 


0.018 


this work (HAO) 


V 


10.521 mag 


0.019 


this work (HAO) 


Ic 


9.788 mag 


0.037 


this work (HAO) 


J 


9.347 mag 


0.023 


Cutri et al. (2003) 


H 


9.069 mag 


0.210 


Cutri et al. (2003) 


Ks 


9.004 mag 


0.020 


Cutri et al. (2003) 


WISEl (3.4 ^m) 


8.932 mag 


0.024 


WISE2 (4.6 ^lu^) 


8.970 mag 


0.022 




WISES (12 /im) 


8.934 mag 


0.030 




WISE4 (22 iira) 


8.662 mag 


0.344 






0.16 


0.04 


this work 


d 


125.1 pc 


4.6 


this work 


Teff 


5879 K 


29 


this work 


log 9 [cgs] 


4.48 


0.15 


this work 


[Fe/H] 


-0.01 


0.05 


this work 


'^micro 


1.40 km 


0.05 


this work 




1.07 Mq 


0.08 


this work 




0.99 Rq 


0.18 


this work 


^systemic 


25.9 km 


0.2 




Vrotsini 


<3 km s-i 







Note. — 



A summary of some of the basic observational properties of TYC 4110-01037-1. 
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Table 3 
Properties of MARVELS-3B 



Parameter 



Value 



Tc (BJDTijB - 2450000) 

P (days) 

e 

u) (radians) 
K (m/s) 
7TiVG (m/s) 
7APO (m/s) 
ecos(a)) 
esin(c<;) 

Tp (BJDtos - 2450000) 
a (sin i = 1) 

M (Mj^p) 



5175.97 ±0.14 
78.994 ± 0.012 
0.1095 ±0.0023 

4.J80_„ (,42 
4199 ± 11 
1338 ± 19 

2945.0l9°g° 

U.UOO(_Q QQ43 

-0.1034 ±0.0027 
5210.321°-^? 
0.38 AU 
>97.7±5.8 



Note. — A summary of some of the basic orbital properties of MARVELS-3B. The quoted minimum mass corresponds to the limiting case 
of sini = 1. Note that Tp corresponds to the time of periastron while Tc corresponds to the time of conjunction. 
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Table 4 

MCMC Properties of MARVELS-3B 



Assumed Prior 


Mass of MARVELS-3B 


Transit Probability 


q (Mass Ratio) 


none (sini = 1) 


97.7 Mj„p 


1 


>0.087 ±0.003 


dN/dq oc q+l 


166.2 Mj„p 


0.0055 


0.149 


dN/dq oc q"^ 


113.0 Mj^p 


0.0129 


0.100 


dN/dq = const 


125.9 Mj„p 


0.0092 


0.112 



Note. — A summary different mass estimates of MARVELS-3B and the resultant mass ratio for different assumed priors on the statistical 
companion mass ratio around G dwarfs and our MCMC analysis. The prior "none" corresponds to the minimum mass of MARVELS-3B, 
simply assuming sini = 1. For all other priors, the quoted companion mass and q values are the median mass determined from our MCMC 
analysis described in Section 4 and shown in Figure 6. 
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Table 5 

Properties of the Candidate Tertiary Companion to TYC 4110-01037-1 



Filter 


Magnitude 


J 


13.59 ± 0.13 


H 


13.04 ± 0.24 


k' 


12.84 ± 0.05 



Note. — 



The computed apparent magnitude of the candidate tertiary companion to TYC 4110 detected in Figure 8 is compiled. 



0.1 



1.0 

X (|im) 



10.0 




0.1 1.0 10.0 

X (jxm) 



Fig. 1. — Top: The observed near UV through IR SED for TYC 4110-01037-1 is shown along with a best-fit NextGen model atmosphere. 
The resultant fundamental stellar parameters from this fit agreed to within 1-cr with the stellar parameters determined from analysis of 
moderate resolution spectra (Section 3.1) Bottom: By constraining Teff ^.nd log (g) to the spectroscopically determined values during the 
SED fit, we are able to better constrain the total line-of-sight extinction to be A^^ = 0.16 ±0.04. Using this Ay, we estimate the distance to 
TYC 4110-01037-1 to be 125.1 ±4.6 pc. 
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6500 6000 5500 5000 4500 

Teff [K] 

Fig. 2. — The observed stellar parameters for TYC 4110-01037-1 (red data) are compared to a Yonsei-Yale stellar (Demarque et al. 2004) 
evolutionary track for a M* = 1.07 Mq star with [Fe/H] = -0.01. Ages of 1.0, 2.0, 5.0, and 6.0 are indicated in blue, and 1-a deviations in 
the evolutionary track arc shown in the shaded region. 
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Fig. 3. — The periodogram for RV measurements of TYC 4110-01037-1 (panel a) indicates the strong presence of a ~79-day periodic 
signature, indicated by the dashed vertical hne. This signal is not seen to coincide with any strong feature in the periodogram for all stars 
located on the same SDSS-III plate as this object (panel b), indicating it is not caused by a known artifact of our instrument or reduction 
pipeline. 
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O 800 1000 1200 1400 1600 



BJD-TDB - 2454000 

Fig. 4. — The derived relative radial velocities from the MARVELS (blue points) and SARG (red points) spectrographs are overlayed with 
the best fit orbital solution described in Section 4 and compiled in Table 3. Residuals to this fit are shown in the bottom panel. 
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6 0.0 0.2 0.4 0.6 0.8 1.0 



Phase 

Fig. 5. — The phase-folded radial velocity curve is shown, for a MARVELS-3B period of 78.994±0.012 days and eccentricity of 0.1095 
±0.0023. Based on the derived stellar mass of 1.07 to°gl Mq (Table 2), we determine the minimum (sin i = I) mass of MARVELS-SB to be 
> 97.7±5.8 Mj„p. 
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Fig. 6. — The cumulative probability that the mass of MARVELS-3B is less than a given mass is shown, for 4 different priors on the 
companion mass ratio described in Section 4. 
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Fig. 7. — The continuum normalized spectrum of HD153458 was subtracted from the continuum normahzed spectrum of TYC 4110-01037-1. 
The l-cr Poissan errors in the data are illustrated with grey boxes. As discussed in Section 4.3, we note that minor mismatches between the 
properties of the reference and science spectra could produce noticable subtraction residuals in line wings, which suggests that only strong, 
repeatable deviations in these difference spectra should be interpreted as real contributions from companion(s). A M3-type companion would 
have contributed a 2% flux enhancement to the system. However, the difference spectrum does not exhibit >3-it deviations above the level of 
the Poisson noise (0.8%) and continuum normalization uncertainties present in the data, which sets the lower mass limit for MARVELS-3B 
that can be ascertained from these specific data. 
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TYC 4110-01037-1 

E 



NIRC2, K' 
January 7, 2012 




Fig. 8. — Keck adaptive optics image of TYC 4110-01037-1 in K . We have detected a faint candidate tertiary companion (indicated by 
the arrow) with red colors that is separated by 986 -I- / — 4 mas from the primary star. If it is physically associated with the primary, it is 
most likely a dM3-dM4 star. 
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Fig. 9. — The mass ratio, q, versus orbital period distribution of binary stars having periods <10,000 years and primary star T^f f < 6000 K is 
shown. Filled blue circles represent K and M dwarf primary stars, and is based on data compiled in literature (Halbwachs et al. 2003; Maceroni 
& Montalban 2004; Burgasser et al. 2007; Becker et al. 2008; Bender & Simon 2008; Rucinski & PribuUa 2008; Dimitrov & Kjurkchieva 2010; 
Davenport et al. 2011) and at http://www.vlmbinaries.org. Amongst these low mass stars, a distinct lack of low q binaries at short orbital 
periods is present. A similar trend is also observed in the distribution of confirmed binaries around solar-like (0.8 Mq < M < 1.1 Mq) stars 
(filled green squares), as seen in data repoduced from tabulated data in literature (Halbwachs et al. 2003; Bender & Simon 2008; Raghavan et 
al. 2010; Sahlmann et al. 2011). The mass ratio of the TYC 4110 (large red cross and an arrow of arbitrary size), OGLE-TR-122 (small black 
cross; adopted from Pont et al. 2005), Kepler-16 (small black cross; adopted from Doyle et al. 2011), and vB 69 (small black cross; adopted 
from Bender & Simon 2008) systems are unique in that they lie within this mass ratio-period deficit. 



